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Abstract 

The asymptotic quasinormal frequencies of the brane-localized (4 + n)- 
dimensional black hole are computed. Since the induced metric on the brane 
is not an exact vacuum solution of the Einstein equation defined on the brane, 
the real parts of the quasinormal frequencies uj do not approach to the well- 
known value Tjj In 3 but approach to Th In k n , where k n is a number dependent 
on the extra dimensions. For the scalar field perturbation Re(uj/T}{) = In 3 is 
reproduced when n = 0. For n/0, however, Re(uj/TH) is smaller than In 3. 
It is shown also that when n > 4, Im(u/Tjj) vanishes in the scalar field pertur- 
bation. For the gravitational perturbation it is shown that Re{ijo/T}{) = In 3 
is reproduced when n = and n = 4. For different n, however, Re(uj/T}{) 
is smaller than In 3. When n = oo, for example, Re(uj/T}{) approaches to 
ln(l + 2cosv / 5tt) ~ 0.906. Unlike the scalar field perturbation Im(uj /Th) 
does not vanish regradless of the number of extra dimensions. 
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The stability problem of the black holes, when perturbed by the external fields, is a 
long-standing issue in the context of the general relativity [1-4]. It is well-known that all 
perturbations are radiated away, which is characterized by the quasinormal modes [5] . These 
quasinormal modes are defined as solutions of the perturbation wave equation, belonging 
to complex-characteristic frequencies and satisfying the boundary conditions for the purely 
outgoing waves at infinity and purely ingoing waves at the horizon, i. e. 

V ~ e iuJZ as z -> -oo (1) 
\p ~ e - iuJZ as z -> oo 

where z is an appropriate "tortoise" coordinate and the time dependence of the fields is 
taken as e luJt . 

In our notation the quasinormal frequencies uj should satisfy Im{uj) > 0. As a conse- 
quence, the quasinormal modes diverge exponentially at both boundaries. This makes it 
extremely difficult to determine the quasinormal frequencies numerically. This is a main 
reason why only very few frequencies with moderate imaginary parts were known [6]. 

About two decades ago Leaver [7] found a possibility to compute the quasinormal frequen- 
cies without having to deal with the corresponding quasinormal modes numerically. Instead 
of the solutions he used the recursion relation which provides an infinite continued fraction. 
However, this method has a technically problem of convergence when Im{uj) » Re{uj). 
This defect of the continued fraction method was mostly removed by Nollert [8] by comput- 
ing the "remaining" infinite continued fraction. 

For the case of the Ad Schwarzschild black hole Nollert showed using his improved numeri- 
cal method that the asymptotic quasinormal frequencies for the scalar field and gravitational 
perturbations become 

uj= U + + 0.0874247 (n = 0, 1, 2, • • •) (2) 
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with, for simplicity, assuming r^ — X where rjj is an horizon radius. This was confirmed by 
Andersson [9] by the phase integral method, which is an improved WKB-type technique. 



After few years Hod [10] claimed surprisingly that the numerical number in Eq.(2) is 
identified as 

0.0874247 -> — = T H In 3 (3) 

47T 

where T# is an Hawking temperature. This identification and the Bohr's correspondence 
principle naturally imply that the minimal quantum area is 4 In 3, one of the values 4 mA; 
suggested in Ref. [11]. This is intriguing from the loop quantum gravity [12] point of view 
because it suggests that the gauge group should be 5*0(3) rather £77(2) [13]. Subsequently 
the identification (3) was analytically shown in Ref. [14,15]. Especially in Ref. [15] the 
authors transformed the boundary condition of the quasinormal modes at the horizon into 
the monodromy in the complex plane of the radial coordinate. We will use this method to 
compute the asymptotic quasinormal frequencies of the brane-localized (4 + n)-dimensional 
Schwarzschild black hole. For the quasinormal modes of the other asymptotically flat 1 black 
holes see Ref. [18-21] and references therein. 

Recently, much attention is paid to the higher-dimensional black holes. Besides its own 
theoretical interest the main motivation of it seems to be the emergence of the TeV scale 
gravity arising in the brane- world scenarios [22-25], which opens the possibility to make a 
tiny black holes factory in the future high-energy colliders such as LHC [26-30]. In this 
reason the absorption and emission problems of the higher-dimensional black holes were 
extensively explored recently [31-34]. 

The lower quasinormal frequencies for the brane-localized 5-dimensional rotating black 
holes were recently computed numerically [35]. In this letter we would like to go further 
for the study of the brane-world black holes by examining the asymptotic quasinormal 
frequencies of the brane-localized Schwarzschild black holes. Since the metric of the brane- 
localized Schwarzschild black hole is induced by the higher-dimensional bulk metric, it is not 
a vacuum solution of the Einstein equation defined on the brane. We will show that this fact 



See also Ref. [16,17] for asymptotically non-flat black holes 
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makes the real part of the asymptotic quasinormal frequencies not to be the well-known value 
T H \n3 but to be T H \nk n , where T H is an Hawking temperature of the higher-dimensional 
black hole and k n is a number dependent on the extra dimensions. 

We start with the (4 + n)-dimensional Schwarzschild black hole whose metric is given by 
[36,37] 



ds% = -h(r)dt 2 + h-\r)dr A + r 2 dtt 2 n+2 



where 



t) 



n+l 



h(r) = 1 - 
dtt 2 n+2 = del + sin 2 Oi 



(4) 



(5) 



d0 2 2 + sin 2 2 {<$! + ••• + sin 2 n (d6 2 n+l + sin 2 6 n+1 d<f) 2 ) • • •} 



and th is an horizon radius. In Eq.(4) the subscript "B" stands for "bulk". For the consid- 
eration of the brane- localized black hole we should consider the Ad- induced metric of ds 2 B , 
which is 



ds 2 4 = -h(r)dt 2 + h~\r)dr 2 + r 2 (d6 2 + sin 2 Gd<j?). 



(6) 



Now, it is easy to show that the equation □$ = 0, which governs the scalar field pertur- 
bation leads to the following radial equation: 



h(r) d 
r 2 dr 



r 2 h(r)—— 
dr 



+ 



u 2 -h{r)- 



+ 1) 



R = 0. 



(7) 



When deriving Eq.(7), we used the factorization condition $ = e wt R{r)Y^^ m {6 , 0). Defining 
R — ^/r, one can show that Eq.(7) reduces to the following Schrodinger-like equation 



^ - V s {r) 



* = 



(8) 



(9) 



where Vs(r) is an effective potential of the scalar field perturbation, which is given by 

Vs(r) = h(r) + 

Now, we would like to derive the wave equation which governs the gravitational pertur- 
bation. In order to derive it we should change the metric itself, say ds 2 — > ds\ = ds\ + 5s 2 
where 



dP 

5s 2 = [H (r)dtd(j) + H^drdQ] e iM sin0-f (cos0) (10) 

do 

with assuming H , Hi « 1 for the linearization 2 . Then, it is straightforward to derive the 
Ricci tensor R^ v and curvature scalar R. 

It is worthwhile noting that since the metric ds\ in Eq.(6) is an induced one from the 
higher- dimensional bulk metric, it is not an exact vacuum solution of 4d Einstein equation. 
Thus, it may satisfy the non- vacuum Einstein equation = where = i? M „ — g tMl/ R/2 
and T^ v is an energy-momentum tensor. Since the perturbation in general changes both the 
geometry of the spacetime and the energy-momentum tensor, adding 5s 2 to ds 2 should 
makes the Einstein equation as £^ v + 5S^ U = T^ v + 5T^ V where 5£^ v and 5T^ V are order of 
H or Hi. Although we can not compute ST^ by conventional method because we do not 
know the matter source nature, we can compute it as following. Firstly, we note that the 
nonzero components of 5T^ U are only ST^, 5T r( p and 5Tg^. This can be easily conjectured 
by computing the Einstein tensor. One constraint is a covariant conservation of the energy- 
momentum tensor (J 7 *" + 5T^ U ).^ = 0. When v = t, r and 9, this is automatically satisfied. 
Thus, this conservation law generates one constraint. Second constraint comes from the 
fact that the effective potential which will be derived later should coincide with the well- 
known Regge- Wheeler potential when n = limit. The final constraint comes from the fact 
that the (t, 0) component of the Einstein equation should be derived from (r, 0) and (9, (j>) 
components. These constraints uniquely determine 5T^ V and the final expressions of the 
Einstein equation are 



1 d if o , 



£(£+!) -2 1 1 dh ld 2 h 

4T 2 + 27 2 W + Yrlh- + Idr 1 



H (11) 



iu (dHi 2 \ 1 h2 (2dh Ad 2 h d 3 h\ ^ 

4 I dr r 1 4iu \r 2 dr r dr 2 dr 3 / 1 



iu , . . / dH 2 TT \ 



+ 1) -2 u 2 7 Al . I dh ld 2 h 
- —h-^r) + — — + - — 



4r 2 4 2r dr A dr 2 



Hi = 



2 The metric change in Eq.(10) corresponds to the odd-parity gravitational perturbation. For 
detail see Ref. [4] 
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, 1 / x „ , / x dH-i dh „ 
iuh^MHo - ^ r — 1 - -r#i = 0. 

dr dr 

Eliminating H in Eq.(ll) appropriately, one can derive a second-order differential equation 
for solely Hi in the form 



+ 3/1 (r) I — 



+ 



dr r J dr 

2N 



(12) 



(f ) i//- 2 (/-)- r v ' 1 ~ + --)/^'-(r) 



7(£+l)-2 { 4dh" 
r dr 



H x = 0. 



Defining Hi = rh 1 (r)^(r), we can derive the following Schrodinger-like equation: 

2 



Hr)^\ ^+[u 2 -V G (r) 



* = 



(13) 



where Vc(r) is an effective potential of the gravitational perturbation, which is given by 

' £(£ + !) Idh 2 



Mr) = Kr) l^r^ + - Tr --^~ Kr)) + ~, 



d^K 
dr 2 



(14) 



Using dh{r)/dr — (n + 1)(1 — h(r))/r, one can express the effective potentials Vs(r) and 
Va(r) as following 



V eS (r) = h(r) 



+ ^(l-h(r)) 



(15) 



where 



n + 1 for scalar field perturbation 

— (n + l) 2 — 2 for gravitational perturbation 

Thus we can treat the scalar and gravitational perturbation in an unified way. 
Now, we define a tortoise coordinate z as 



(16) 



d _ M d 



1 - 



n+1' 



d_ 

dr 



(17) 



Integrating Eq.(17), we can make an explicit expression of z in the form 



z = r + 



n + 1 



J2 e i27Tj/(n+l) j n 



3=0 



(18) 



When deriving Eq.(18), we fixed the integration constant by imposing that r = corresponds 
to z — 0, which is crucial [15] for the calculation of the asymptotic quasinormal frequencies. 
From this tortoise coordinate and the boundary condition (1) it is straightforward to compute 
the monodromy around the singularity r = 1 

M(l) = e u/TH (19) 

where T# is an Hawking temperature of the (4 + n)-dimensional black hole, i.e. Th = 
(n+ l)/47rr H . Following Ref. [15] we will compute another expression of the monodromy by 
following a contour given in Fig. 3 of Ref. [15] and using an analytic continuation. Equating 
these two expressions enables us to calculate the asymptotic quasinormal frequencies. For 
the following calculation it is important to note that near the naked singularity r ~ the 
tortoise coordinate z is proportional to r n+2 as 

s n+2 

n + 2 \r H/ 



Adopting a method used in Ref. [15], we can now compute the asymptotic quasinormal 
frequencies. The explicit expression of the effective potential given in Eq.(15) is 



1- ( — 

r 



r H \n+i 



+ 1) + ^(rH\ n+1 



(21) 



Near the naked singularity r = 0, therefore, the effective potential can be written by its 
dominant term 

Thus, the wave equation 

d 2 ^ 



dz 2 

simply provides a solution 



+ P " ^eff(O) * = (23) 



^o :A ~ V2mJz [A + J v (uz) + A_J_ v {ujz)\ (24) 

in this regime where 



V4 (n + 2) 2 ' 

Taking a ujz — > oo limit in ^o,A, we can obtain one asymptotic solution 



(25) 



oo.A 



^ £ in(l+2v)/4 + A _ e m(l-2u)/4 



(26) 



with a constraint 



^ -mt(1+2i/)/4 + A _ e -iir(l-2v)/4 = q 



(27) 



To follow the contour given in Ref. [15] we should turn to an angle 3/r around z = in 
\I/o,a- Then we can derive another solution in the r ~ regime 



^o,b ~ V27uJze^ A + J u (ujze i37r ) + A^J_ u (ujze 



(28) 



Since J v (z) = z u ip(z) where ip(z) is an even analytic function in the entire region of the 
complex plane, it is easy to show that \l/o,s reduces to 



^o,b ~ e i3w/2 V2nujz 



A + e i3ni/ J u {uz) + A_e- l '^ u J_ v {ujz) 



—iZ-Kv 



(29) 



Taking ujz — > — oo limit in ^o,b, we can obtain another asymptotic solution 



oo,B 



A+( J^(l+2u)/4 + ^^ e «57r(l-2z/)/4 g-iwz _|_ ^ +e i7ir(l+2i/)/4 _|_ ^_ e i77r(l-2i/)/4 



Now, we can follow the contour to return to the initial point through the asymptotic 
region. Since the second term in Eq.(30) is exponential small in this region, the monodromy 
can be written as 

A i^{l+2v)/4 + A _ e i5n(l-2u)/4 



M ^ ^ +e m(l+2i/)/4 + A _ e i-K{l-2u)/4 ■ 

With an aid of the constraint (27) this simply reduces to 

M(l) = -(1 + 2cos2ttz/). 



(31) 



(32) 



Thus equating Eq.(19) with Eq.(32) gives 



e u/T H = _(1 +2COS27TZ/). (33) 

Now, we turn to the scalar field perturbation. In this case v becomes 

V8 = )lz--fr+W- (34) 

Note that v s = 2/3 when n = 4, which makes the rhs of Eq.(33) to be zero. When n < 4, 
the asymptotic quasinormal frequencies become 

^ = 2ni(n + l-) + ln(l + 2 cos2ttz/) (n = 0, 1, 2, • • •). (35) 

Tff V 2/ 

Of course, v s = when there is no extra dimension, which gives a well-known result 
Re^us/Tn) = In 3. When extra dimensions exist, the real parts of the asymptotic quasinor- 
mal frequencies are summarized in Table I 3 . 



n 


Re(u s /T H ) 





ln3 = 1.09861 


1 


ln2 = 0.69315 


2 





3 


ln((3 - V5)/2) = -0.96242 


4 


— oo 


Table I 



It is interesting to note that Re{ujjTu) vanishes when there are two extra dimensions. 
When n > 4, the rhs of Eq.(33) becomes positive. Thus the asymptotic quasinormal fre- 
quencies become 



3 In Ref. [35] it is shown that when n = 1, the real part of the lowest quasinormal frequency is 
0.273 which is much larger than the asymptotic value X#ln2 = 0.055 
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ln(— 1 — 2 cos27rz/5) 



(36) 



In this case, therefore, we get Im{ujs /Th) = 0. If there are infinite extra dimensions, vg 
becomes 1/2, which makes Rc{ujs/Th) = 0. 
In the gravitational perturbation v becomes 



It is interesting to note that v G — 1 when n = and n = 4. Thus in addition to the case of 
no extra dimension we get a well-known result Re(u>G/T H ) = In 3 when there are four extra 
dimensions. Unlike the case of the scalar field perturbation, however, the rhs of Eq.(33) is 
always negative regardless of n. Thus, the asymptotic quasinormal frequencies are always 
given by 




(37) 





(38) 



The real part of ujg/Th is plotted in Fig. 1. 
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FIGURES 




In 3 



FIG. 1. The n-dependence of Re(ujg/Th)- It is interesting to note that the well-known result 
Rs(ujg/Th) = hi3 is reproduced when n = and n = 4. For other value of n Re^c/Tn) is shown 
to be smaller than In 3. 



Fig. 1 shows that ReioJc/Tn) becomes the well-known result In 3 when n = and 
n = 4 as mentioned before. When n is other number 4 , Re(u G /T H ) becomes smaller than 
In 3. In particular, if there are infinite extra dimensions, v G becomes y/5/2, which gives 
Re(u G /T H ) = ln(l + 2 cos VEtt) « 0.906. 

In this letter we have computed the asymptotic quasinormal frequencies of the brane- 
localized (4 + n)-dimensional Schwarzschild black hole. For the case of the scalar field 
perturbation the well-known result Re(u>s/T H ) = In 3 is reproduced when there is no ex- 
tra dimension. When < n < 4, Re(us/T H ) is shown to be smaller than In 3 When 
4 < n, Im(us /Th) becomes zero. In the case of the gravitational perturbation, however, 

4 When n = 1, Ref. [35] shows that the real part of the lowest quasinormal frequency is 0.805 
which is much larger than the asymptotoc value 0.086. 
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Itti^g/Th) does not vanish regardless of the number of the extra dimensions. In particular, 
the well-known result Re{uJa/T H ) = In 3 is obtained when n = and n — A. 

It is of interest to extend our computation to the brane-localized (4 + n)-dimensional 
Reissner-Nordstrom and Kerr black holes. It is known that the asymptotic quasinormal 
frequencies of the Ad non-extremal Reissner-Nordstrom black hole have real part T H In 2 
unlike the Schwarzschild black hole case [15]. It is of interest to examine how the existence of 
the extra dimensions changes the quasinormal frequencies. In the case of the Kerr black hole 
the superradiance effect seems to play an important role [19]. Since it is proven generally that 
the superradiance modes exist in the brane-localized and bulk rotating black holes [38,39], 
it seems to be interesting to examine how these modes affect the asymptotic quasinormal 
spectrum. It is also interesting to confirm our result by adopting an appropriate numerical 
technique. 

The most interesting issue at least for me is to explore how the n-dependent quasinormal 
modes affect the Hawking radition of the black holes. Similar issue was studied long ago by 
York [40]. We would like to study this issue in the future. 

Acknowledgement: This work was supported by the Korea Research Foundation under 
Grant (KRF-2003-015-C00109). 
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